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Abstract: We compute a T-odd triple vector correlation for the /r —)• eq decay and the 
/X —>■ e conversion process, hnding simple results in terms of the CP violating phases of the 
effective Hamiltonians. Then we focus on the minimal Left-Right symmetric extension of 
the Standard Model, which can lead to an appreciable correlation. We show that under 
rather general assumptions, this correlation can be used to discriminate between Parity or 
Charge-conjugation as the discrete Left-Right symmetry. 
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1 Introduction 

Lepton Number Violating (LNV) and Lepton Flavor violating (LEV) processes are forbid¬ 
den in the Standard Model (SM) and are thus a good probe of new physics. In principle new 
physics brings also new sources of CP violation and therefore time reversal (T) symmetry 
violation in any local, Lorentz invariant quantum field theory. 

Motivated by this we explicitly compute T-odd triple vector correlations for the LEV 
/X — 7 - 67 decay and /x —)■ 6 conversion process, since much of the present and future experi¬ 
mental efforts are devoted to these two processes. The MEG collaboration reports the best 
experimental limit for the /x —)• 67 decay [ 1 ] 

Br(/x ^ 67 ) = < 5.7 X 10-13 ( 1 . 1 ) 

r(/x eUf.l'e) 

and the SINDRUM II collaboration gives the strongest limits for the /U —>■ 6 conversion 
process [2, 3], namely 

Br(;u + Ti(Au)) ^ 6 + Ti(Au)) = ^ < 6.1(7) x lO'^^, ( 1 . 2 ) 

-L capt 
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where Fcapt is the muon capture rate in the vicinity of a nucleus. Upgrades of ongoing 
experiments have been considered with the hnal goal of achieving a sensitivity around 
10“18 — 10“^® [4-7]. Given the current limits and the future improvements, there exist the 
possibility of having enough statistics to start probing CP violation beyond the SM in the 
next round of experiments. This is suggested and studied in [8, 9]. 

In this work we focus on quantities that test T violation in the absence of hnal-state 
interactions and among these quantities are triple vector correlations made up of the mo¬ 
menta or spins of the participating particles [10]. In [11], it is suggested that triplet vector 
correlations can be used to probe CP violation in the /x —)> e conversion process. Here we 
present the hrst analytical computation for the correlation suggested in [11] for the n ^ e 
conversion process and we extend their work in two ways: hrst, we compute the correlation 
for the /X —>■ ey decay and second we include the full set of effective operators that enter 
the /X —>■ e conversion process. 

In section 2 we introduce some theoretical tools for the /x —)• ey and the /x —)• e conversion 
process. In section 3 as an example of a theory that gives order one contribution to the triple 
vector correlation, we briehy introduce the minimal Left-Right (LR) symmetric extension 
of the SM. In section 4 and 5 we present the analytical computation of the triple correlation 
in the /x —?• eq and /x —e conversion process respectively. Then in section 6 and for both 
processes, we study these correlations in the context of the minimal LR model, for both 
parity and charge-conjugation as the LR symmetries. Finally in section 7 we present our 
conclusions. 

2 General theory 

2.1 /X —eq process 

The /X —)• eq decay is predicted to be negligible small in the SM with massive neutrinos, 
therefore if this process is seen it implies that new physics is behind it. The effective 
Hamiltonian for this process is of the form 

= ^^^^e{pe)a^uF>^^ALPL + ArPr)p{p^) + h.c., ( 2 . 1 ) 

where e is the electromagnetic coupling constant, is the electromagnetic held strength 
for the photon held, Gr \s the Fermi constant, P[r^l) = 2(1 4= qs) , is the muon mass 
and e{pe) and p{p^) are the spinors for the electron and muon respectively. For this process 
we use the gamma matrices in the Weyl basis and the coefficients Ar and Ar are calculated 
within a given physical model. 

2.2 p ^ e conversion. Theory and Effective Hamiltonian 

Theoretical studies of this process were performed in the past [12-15]. In [15] the outgoing 
electron coming from the conversion process, belongs to one of the states in the continuum 
energy spectrum for the Coulomb potential and as a matter of fact the outgoing electron 
must be treated as a plane wave. One way to argue this is by noticing that an electron in 
the continuum energy spectrum, is described by a Dirac spinor in the angular momentum 
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basis. Experimentally, the detected electron has a define 4-momentum implying that the 
outgoing electron must be a plane wave. 

In this work we present a method for computing a triple vector correlation that tests T- 
violation in the /r —)• e conversion process for various nuclei. We make use of the formalism 
developed in [16]. 

We use the following representation for the 7 matrices 


and 


70 = /3 




(^ 111 / — 


1 

2 




75 = -*71727370, 


( 2 . 2 ) 


(2.3) 


where the cjj are the Pauli matrices where i = 1,2,3 and the index ^ takes the values 
// = 0,1,2,3. 

The Dirac’s equation for the central field problem in polar coordinates is given by (the 
energy is given in units of the electron mass) 


E^Ij = H^Ij = [-i75S.(|- + --^K) + V + /3]^, (2.4) 

or r r 

where 

^ = ^[7i,7fc] cyclic). (2.5) 

i 

K = L + l). (2.6) 


V is the Coulomb potential and L is the orbital angular momentum. 
We write the wave function as [17] 




( 9K{r)XK \ 


(2.7) 


such that K^|:^ = —ki/’k and = /Ui/’k, where J 3 is the third component of the total 

angular momentum J. The radial functions and obey the differential equations 


= -—9.{r) + {E-V + l)Ur), (2.8) 

dr r 

= —Ur) -{E-V- l)gU)- (2.9) 

dr r 

In the high energy limit -all the masses are set to zero- and from eqs.(2.8) and (2.9), 
fU) and gU) satisfy 

, 9—k — fn- 


3 - 


f—K - 9k 


( 2 . 10 ) 






From here on we make use of this result for the spinor '4>^e describing the electrons 
coming from the conversion process. The initial muon instead is described by with the 
quantum numbers, /r = and k = — 1 and we choose the normalization 

J = 1 . ( 2 - 11 ) 


For the electrons in the continuum-energy states we use the same normalization con¬ 
sidered in [15], namely 


j - e'). 

In the conversion process the effective Hamiltonian is given by [15] 


( 2 . 12 ) 


Heff = PieE^ty + PneE^t, + h.c.) 

Gp 

+ ^ X] \.^9LS(q)^PRd‘ + gRS{q)^PLd)m + {9LP{q)ePRd‘ + 9 RP{q)^PLd)qi5q 

^ q=u,d,s 

{9LV{q)ei^PL9 + 9RV{q)e'y^PR9)qit^q + {9LA{q)n^PL9 + 9RAiq)ej^PR9)qit^l5q+ 
\{9LT(q)e(y^'' Pr 9 + 9RT(q)e(T^''PL9)qcr^,uq + h.c.]. (2.13) 


The nuclear form factors were calculated in [18]. The wave function for the muon and the 
electrons in the presence of a central held were obtained in [14, 15]. In particular in [15] 
updated data for the proton and neutron densities were used. 

In the limit of r —?• oo it can be shown that the general solution for a Dirac particle in 
a Coulomb held at hrst order in Ffe// is of the form [16] 


V’as 





^/E + lx'iip) 


+ (2.14) 


where p is in the direction of the outgoing electron. The phases are the usual ones 
appearing in scattering problems in the presence of a Coulomb held and are given by 

6k = y\n2pr - argr (7 -k iy) + Pk - ^vry, (2.15) 

y = aZE Ip, 7 = k? — a^Z'^, =-—— (2.16) 

j + iy 

where Z is the atomic number, a = e^/dvr and p is the modulus of the 3-momentum p. 
We consider states with k = ±1, hence the only term relevant for our discussion is -the 
remaining ones are just an overall phase in the solution tpas- 
Finally the total conversion rate per unit hux is ^ 

(^conv = P f dn-tpl^'lljas = ^ ^ | (V'^|FFe//|' 0 i) P- (2.17) 

^See appendix B for a more detailed discussion on this issne. 
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3 The minimal Left-Right symmetric theory 


As an example of a complete and predictive theory of lepton number violating phenomena 
we consider the minimal LR symmetric extension of the SM [19-22], In this model the 
gauge group is SU{2)l x SU{2)pi x U{1)b-l with an additional discrete symmetry that 
may be generalized parity (P) or charge conjugation (C) -for reviews see [23-25], It relates 
the smallness of neutrino masses to the near maximality of parity violation al low energies 
through the seesaw mechanism [27-31], The scalar sector contains the following fields [26- 
29] 


$ = 




^L,R 


( ^L,r/ ^ 

[ sin J ’ 


(3.1) 


where is in the (2,2,0) representation of SU{2)l x SU{2)pi x U{1)b-l and the two scalar 
triplets Al and Ab, belong to the (3,1,2) and the (1,3,2) representations respectively. The 
Yukawa interactions of leptons with the scalar triplets have the form 


Cy = Tl(Tce.$ + Y$I>)Lr + UlICW^^AlLl 

AjjL/j) + h.c., (3-2) 


$ = (T 2 <h*(T 2 , Ll is the lepton doublet of the standard model (L|( = {y 1 )l) and Lr is its 
right-handed analogue that we denote as = {N l)ji where N is the heavy Majorana 
neutrino. The Ya is the Yukawa coupling of the field a, where a = {$, A^j}, 

Under the discrete left-right symmetry the fields of the theory transform as: 


P : 


Pflp-^ = lofR 


PAl,rV ^ = -Arl 


C : 


C/lC-^ = CUnf 


(3,3) 


where the usual charge conjugation operator is given by C = Z 727 o- 

Invariance of the Lagrangian under the LR symmetry leads to the following relations 
between the Yukawa couplings of the theory, namely 


P : 


^^R,L ^^L,R 

= yI 

Y.s> = Y^ 


C : 


Ya 


R,L 


= YX 


Y^ = Y^ 

Ke, = Y^ 


(3,4) 


In the mass eigenstate basis the flavor changing charged current Lagrangian is given 
by 

Ccc = -j={vLVlM^ lIl + NrVI\^rIr) + h.c., (3,5) 

Vr is the right-handed analogue of the PMNS mixing matrix Vr- In general it has three 
different mixing angles and six arbitrary complex phases and we parametrize it as 


Vr = K.VrKm, 
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with Ke = diag(e*'^=, Kp^ = diag(l, The matrix Vr has three mixing 

angles and the dirac phase 5. We choose for Vr the standard form for the CKM matrix 
shown in the PDG [32], 

The interaction terms of charged leptons with the doubly-charged scalars are 

£a = + hlcYij++ii, + h.c., ( 3 . 6 ) 

If charge conjugation is the discrete LR symmetry, the charged lepton masses are 
symmetric and in this case the Yukawa couplings in (3.6) satisfy (for reviews on this topic 
see references [23-25]) 

= (YkJ- (3.8) 

For parity and in the more interesting phenomenological situations, the charged lepton 
masses matrices are almost hermitian [33]. In [34] it was realized that it implies the near 
equality between the Yukawa couplings shown in Eq. (3.6) i.e. 

= Y^^ + O(tan 2/3 sin a). (3-9) 

The vacuum expectation values of the neutral fields belonging to <I> are such that = vi 
and {(P 2 ) = '^ 26 *“) where /3 is the ratio v^jvx and a is the spontaneous phase. In [35, 36] 
it is shown that tan 2/3 sin a < 2rrn,/mt {rrih and are the bottom and top quark masses 
respectively), so that the Yukawa coupling of the doubly charged scalar are nearly equal 
[34]. 

It is a remarkable feature of the minimal LR theory, that the TeV energy scale accessible 
at the LHC through the Keung-Senjanovic (KS) process [37] -and its associated LNV 
and LFV, predicts the rate for the neutrino-less double beta decay and low energy LFV. 
This deep connection and the related phenomenology are illustrated in [38, 39]. All these 
processes depend in a crucial way on the elements of the leptonic right-handed mixing 
matrix Vr, for which all its mixing angles, the Dirac phase and two Major ana phases can 
be determined at the LHC [34]. Useful information can also be obtained from EDM of 
the neutron and such [40-44]. This is deeply connected to the study of the strong CP 
parameter, which in the mLRSM turns out to be calculable [45-47]. 

Recently the CMS collaboration [48] has reported an excess in the ee-channel for this 
process at 2.8a, but they claimed that this excess cannot be accommodated in the minimal 
version of the theory -assuming diagonal mixing in the right-handed leptonic sector and 
degenerate masses for the heavy neutrinos. Several works have been proposed [49-52] 
in order to explain this excess and the conclusion was that it would need a higher Left- 
Right symmetry breaking scale, or a more general mixing scenario with pseudo-Dirac heavy 
neutrinos. 

4 Computation of a triple vector correlation in the /r —eq decay 

T-odd asymmetries in the // —)• eq were considered in the past. In [8, 9], it was shown that 
by studying the polarization of electron and the photon coming from the muon decay it is 
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possible to extract the CP-violating phases from the experiment. The conclusion was that 
in order to extract the CP-violating phases both electron and photon polarizations must 
be measured. In this paper instead, we present an alternative way of extracting the CP- 
violating phases of the effective Hamiltonian in the /r —)• ey decay. This is complementary 
to the work presented in [8, 9]. The novelty is that no measurements of the final photon 
polarizations are needed. We consider the T-violating triple vector product 


s^+ ■ (pg+ X Sg+) = cos d* sin 0s, 


(4.1) 


where 9s is the angle between the polarization’s direction (sg+) of the positron and its 
momentum’s direction Pg+, is the angle formed between s^+ and the direction defined by 
Pe+ ^ ^e+ is the azimuthal angle. In Fig.l the reference frame and setup are shown. 

Notice that this quantity changes sign under parity and naive time-reversal transformation T 
defined by t —>■ —t. For processes whose interactions are characterized by a small coupling, 
it can be shown at first order that the connected part of the S-matrix is hermitian [10] 
and therefore the violation of the T symmetry amounts the violation of the time-reversal 
symmetry. 

We define the triple vector correlation as 


(Su+ ■ (Pe+ ^ 'Se+))‘I’ — 


N{cos ^ > 0) — N (cos ‘h < 0) 


^total 


d4»dr/d4» • sgn(s^+ • (pg x Sg+)) 

^total 


(4.2) 


where Ftotal and Ntotal are the total decay rate and the total number of events for the 
initially polarized muon respectively, N (cos <I> > 0) and N (cos <I> < 0) are the number of 
events satisfying cos ‘h > 0 and cos $ < 0 respectively. 

The 4-momenta of the participating particles in the rest frame of the muon are given 
by 


= (m^, 0,0,0), (4.3) 

Pe+ = {Ee, \Pe+\sin9s, |Pg+ | COS 0^, 0), (4.4) 

P 7 = (-E^7>“IPe+|sin0s,-|Pe+|cos0s,O) (4.5) 

where the mass of the positron has been neglected. The energy £'g+ of the positron and 
the energy of the photon are given by 

77 ? 

Ea+^E, = \po+\ = ^. (4.6) 

From the effective Hamiltonian in eqn. (2.1) and eqns. (A.l), (A.4) and (A. 5) in 
appendix A, a straightforward computation gives the following value for the correlation 

/- - \\ ■ /I 3‘"i(AlAp) 

(■5;,+• (Pe+ X Sg+))<!, = sm0*|-^-|2-q-^. (4.7) 


The main advantage of this quantity is that no measurements of the photon polarizations 
are needed. 





y 



Figure 1. Reference frame and the setup for the /r —)■ 07 decay. 


In summary we find that given a source of polarized anti-muons, by measuring the 
3-momentum Pg+ of the outgoing positron and its polarization Sf .+, the asymmetry shown 
in eqn. (4.7) is sensitive to the CP-violating phases of the effective Hamiltonian shown in 
(2.1). In [53-57] it is shown that measurements of the polarization of electrons coming from 
the muon decay are feasible. We assume a 100 % polarized muon flux so that our results 
must be trivially rescaled by the actual polarization of the initial muons. 


5 Computation of a triple vector correlation in the /r — e conversion 
process 


Following the same lines of the last section, we define an asymmetry given by comparing 
the number of events with • (pf, x Se) > 0 with the ones satisfying ■ {pe x Sg) < 0 in 
the /X —>■ e conversion process and we define it as 


{s^l ■ {Pe X Se))$ = 


N{cos ^ > 0) — N (cos < 0) 

^total 

UJconvicOS 4> > 0) - UJconvicOS 4> < 0) 
^r.nn.u 


(5.1) 


where ojconv is the total conversion rate and as previously, is the angle between the plane 
formed by the vectors Pe and Sg and the polarization of the muon We used the same 
coordinate system shown in Fig.l but clearly there is no photon coming from the muon 
decay. 
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A direct computation of the asymmetry shown in Eq. (5.1) gives ^ 


(% • {pe X 


1 . , 

2®“^|Cl|2 + |Ck|2 


+ 0{aZ) + 0{^). 


rUe is the electron mass and 


Cr = DAr + + ~9lv) + + th 

Cr = DAr + + ~gtl) + 


(5.2) 


(5.3) 

(5.4) 


where 

9rS,RS = X/ 9LS,RS{q)^ 9rs,RS = X] (5-5) 

9 9 

9rV,RV = ‘^9lV,RV{u) + 9LV,RV{d)-, 9rV,RV = 9lV,RV{u) + ‘^9LV,RV{d)- (5-6) 

D, are nuclear constants already calculated and tabulated in [15] for various 

elements. and are obtained from the scalar matrix element [15, 18] 

(iV|gg|Af) = + iA- (5.7) 


Z and A are the atomic and mass number respectively, and are the neutron and 
proton densities inside the nucleus. The expression obtained is valid for non-relativistic 
muons and we droped terms of the order aZ and me/ife. In practice equation (5.2) must 
be multiplied by the polarization of the initial muons, which is of the order of 15% in the 
conversion process [58]. 


6 Triplet vector correlation in the minimal Left-Right theory 

As a concrete example of a theory beyond the SM that gives order one values for the T-odd 
triple vector correlation [11] we consider the minimal LR symmetric model. In what follows 
we analyze separately the contributions to the asymmetries (4.7) and (5.2) in the case of 
P and G as the LR symmetries. In [11] it is found that this contribution can be of order 
one, since there are new contributions coming from interactions of charged leptons with the 
singly-charged and doubly-charged scalar fields. 

6.1 p e'j decay 

In this section and for the /U —)• eq decay, we study the contributions to the triple vector 
correlation for both Parity and Charge Conjugation as the LR symmetry. 

Parity as the LR symmetry: in [59] the authors presented a complete study of the 
contributions to several LEV processes in the context of the minimal LR extension of the 
SM and it is found that the branching ratio for this process is of the form 

Br(/i eq) = e^{\AR\^ + |ArP) 

^For more details see section B.2. 
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Figure 2. Plot of the loop function 5*3(a;). 


where 




Al 


Xr, 


167r2 Z^(^Rhn[VRM^2 ^3[Xn) 3 ^2 

n 




1 _ 1 /\/f2 1 /\/f2 

'^i^R)en{VR)niJ.Xn[--j^^ - ^ ] + C>(tan 2/3 sln«), 


-5^ 


( 


Mjv 


53(x) = -- 


X 1 + 2x 


+ 


2>x^ 


X 


8 (1 — x)2 4(1 — x)2 ^(1 — x)2 


(1 - X + logx) + 1]. 


( 6 . 2 ) 

(6.3) 

(6.4) 


Mat^ are the heavy neutrino masses where n = 1, 2, 3. Mw is the W boson mass, is 

the Wr boson mass, Mrr+ is the mass of the heavy scalar Hf and Mc++ are the masses 

(L,R) 

for the left and right doubly charged scalars respectively and finally we use Mi, to denote 
the light neutrino masses. 

Notice that the loop function S 3 is always small as far as M^r is not much bigger than 
Mw^, so that the term with the loop function can neglected for a wide range of the heavy 
neutrino masses (see figure 2) and therefore the correlation defined in (4.7) is suppressed. 
Finally we neglect the contribution of the charged Higgs since its mass cannot be lower 
than (15-20) TeV [35, 60]. This poses no problem for the theory, since its mass emerges at 
the large scale of symmetry breaking [26, 61]. The gauge boson and doubly-charged scalar 
masses can be obtained at the LHC through the so called KS process and the decays of the 
doubly charged scalars [37] in addition with all the mixing angles and the Dirac phase in 
Vr [34]. This is an example of the complementary role played by the high and low energy 
experiments in the establishment of the LR theory [38, 62-67]. 

For the sake of illustration, imagine that type H see-saw is the dominant source of 
neutrino masses i.e. ^ ~ ^R- * 2 ase it is possible to show that the 
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, rr,. 
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Tj 
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0.5 
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NH 

0.0 


IH 


“w„ 




Figure 3. Plot obtained by considering the MEG bound shown in Eq. (1.1). (Right) Normal 
hierarchy case (NH). (Left) Inverse hierarchy case (IH). The colored region is the allowed one. 
(Top) Mass of the heaviest right-handed neutrino =0.1 TeV. (Bottom) Mass of the heaviest 
right-handed neutrino = 1 TeV. 


heavy neutrino masses satisfy the relation [38] 








Ml - Ml 


^2 






- 


~ zbo.03, 


(6.5) 


where the ± corresponds to normal (NH) and inverted (IH) neutrino mass hierarchy re¬ 
spectively. In what follows we denote Mtvq the lightest right-handed neutrino mass, 
the heaviest right-handed neutrino mass and 6 is the Dirac phase present in Vr. In Fig. 3 
and for the two representative values of M^^^ =0.1 TeV and Mri^j = 1 TeV we show the 
allowed region obtained from the MEG bound in the (^Dirac} plane, for both normal 

and inverted neutrino mass spectrum. The region between these values gives rise to the 
exciting LNV signals at the LHC trough the KS process. We assume = 3.5 TeV 

and common masses for the doubly charged scalars M^++ = M^++ = M^ = 1 TeV. The 
reader may ask about the very different behavior obtained for the two values of the heaviest 
neutrino mass chosen, and the point is that this can be readily understood by noticing that 
the amplitude is approximately proportional to [AM^gl = so that a bound 

is obtained for lAM^gl rather on the lightest neutrino mass itself. 

In figure 4 (top) we plot the absolute value for the triple vector correlation given in 
(4.7) in the (Matq, <5)-plane, where one may see that the values of the correlation (4.7) goes 
from 10“® to 10“® in the allowed region. 
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Figure 4. (Top) Contour plots illustrating the absolute value of the asymmetry defined in (4.7) 
as a function of the lightest neutrino mass and the Dirac phase <5 for P as the LR symmetry. 
(Bottom) Contour plots illustrating the value of the asymmetry defined in (4.7) as a function of 
the lightest neutrino mass MjVp and the Dirac phase S (assuming = 0) for C as the LR 

symmetry. (Left) Normal hierarchy for neutrino masses. (Right) Inverse hierarchy for neutrino 
masses. We take the gauge boson mass = 3.5TeV, the heaviest right-handed neutrino mass 

= ITeV and common masses for the doubly charged scalars of Ms = 1 TeV. The mixing 
angles are 6 »i 2 ~ 33.6°, 623 ~ 41.9°, 6*13 ~ 8.7°. 


One would be tempted to conclude that the triple vector correlation may be bigger for 
general values of neutrino masses and mixings. However from eqns. (3.9), the contribution 
to the triple vector correlation shown in (4.7) is bounded to be less 10“^ since tan 2/3 sin a < 
10“^ from the quark masses [35, 36, 47]. The point is that for charged leptons masses (Mi) 
bigger or equal than the Dirac mass of neutrinos (M^), the mass matrix of the charged 
leptons is nearly hermitian leading therefore to nearly equal leptonic left and right mixing 
matrices. This is in complete analogy to the situation in the quark sector studied in [35, 40]. 
Of course it is possible to assume that Md > Mi, but we will not pursue this possibility 
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since in this case the original see-saw mechanism would lose its meaning and one would 
have to invoke accidental cancellations in order to explain the smallness neutrino masses. 
Charge conjugation as the LR symmetry: from eqn. (3.8) we have that 




Al 


107^2 *S'3(X„) 

n ^Wr 

Y.^yR)eniVR)nRXn[--jY^ 
n 5++ 


Xn 

°R 

1 Vir, 

24 m2 , ' 


( 6 . 6 ) 

(6.7) 


Notice that some of the external phases appearing in Vr do not cancel in (4.7) and the 
triple vector correlation is proportional to so that the triple vector correlation is 

not suppressed by the small ^13 mixing-angle. In Fig. 4 (bottom) we show the absolute value 
of the triple vector correlation in the (M^rp, (I)-plane. We take — (pe) = 0 in both normal 
and inverted neutrino mass hierarchies. For {(p^ — (pg) = 7r/4 it will reach in maximum value 
of around 0.5 in almost all the parameter space 

Finally from Fig. 4 (bottom) we conclude that C as the LR symmetry gives larger 
contributions to the triple vector correlation and this because in the parity case, the triple 
vector correlation is suppressed due to the near equality between the Yukawa couplings. 

The bottom line is that in the most interesting region of the parameter space, a value 
for the triple vector correlation bigger than 10“^ can only be the consequence of C as the 
LR symmetry. 

One may ask whether this value for the asymmetry of could be measured in forthcoming 
experiments. Suppose that /i —?• eq is found to be of the order of 10“^'^. In the best scenario 
due to the future experimental improvements on the sensitivity, it would become possible 
to observed at most 10^ events and out of these events one has to select the ones that have 
0s 7^ 0 or 05 7^ TT. Moreover suppose that only the events satisfying 7r/6 < 0^ < 7r/3 may be 
identify in the experiment due to its intrinsic sensitivity. This would imply that we end up 
having sinOgdOg ~ 10^ events in the most optimistic situation. Hence this naive 

argument allow us to conclude that in most optimistic scenario, an asymmetry of the order 
10“^ or bigger would probably be seen in the next round of // —)• eq decay experiments. 


6.2 ^ ^ e conversion process 

In this section we consider the triple vector correlation for the fi ^ e conversion process in 
the context of the minimal LR symmetric extension of the SM where the relevant branching 
ratio is given by [59] 


Br(|U —^ e) 


2(^2^y(p)2 

Leapt 




( 6 . 8 ) 


The values of the capture rate Leapt are tabulated in [68] for several elements. In [59] 
it was shown that the contribution of the doubly-charged scalar may dominate due to a 
logarithmic enhancement and in this case the functions and may be written as 

^(L,R) - 1287r2H(i^7j) log(m2). (6.9) 
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Figure 5. Plot obtained by considering the SINDRUM II bound for Titanium shown in Eq. (1.2). 
(Right) Normal hierarchy case (NH). (Left) Inverse hierarchy case (IH). The colored region is the 
allowed one. (Top) Mass of the heaviest right-handed neutrino =0.1 TeV. (Bottom) Mass of 
the heaviest right-handed neutrino = 1 TeV. 


For completeness we show in Fig. 5 the allowed region obtained by considering the SIN- 
DRUM bound for Titanium shown in Eq. (1.2) assuming the same values for the heavy 
neutrino masses of the last section. From Eq. (6.9) and assuming that the dominant terms 
are the logarithmic enhance ones, the amplitude for the conversion process and the /r —)• eq 
decay are proportional. Therefore a similar qualitative behavior is obtained. We can see 
that the bound obtained is similar to the one of the /U —>■ eq experiment and this is due to 
the fact that the logarithmic enhancement in Eq. (6.9) compensate the a suppression in 
the conversion rate [59]. For Gold the bound one would obtained is similar since the ratio 
between the conversion rates for the two elements is around 0.83. On the other hand, for 
the gold atom relativistic effects of the muon becomes relevant, so that the result shown in 
Eq. (5.2) cannot be trusted in this case. 

Finally the asymmetry defined in Eq. (5.2) takes the form 


_ sinOs _ sing, ^m{ALA\) 

2 \pi'r )\2 2 \Al\‘^ + \Ar\‘^ 

' 1j ' ' li ' 


( 6 . 10 ) 


where it can be seen that this asymmetry has the same flavor structure of the coefficients 
Al and Afi defined previously for the ^ eq decay, therefore the same conclusion obtained 
in the ^ eq case holds for the // —)• e conversion process as well. 
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Regarding the expected sensitivity for the conversion process the arguments we used 
in the /r —)• ey decay apply, but with the difference that the final sensitivity is rescaled by 
a factor of the order of 10“^ due to the depolarization -around 15%- of the muons in the 
conversion process [58]. 

7 Conclusions 

We derived analytical expressions for a T-odd triple vector correlation in the /r —)• ey decay 
and the /x —)> e conversion process and found simple results in terms of the CP-violating 
phases of the effective Hamiltonians. The expression obtained in the fx ^ e conversion 
omits relativistic corrections for the muons, but is otherwise complete. For the /x —)• ey 
decay we conclude that in order to extract the CP violating phases of the theory from the 
experiment, no measurements of the photon polarizations are needed. 

Then as an example of a theory that leads order one values for the triple vector cor¬ 
relation we consider the TeV scale, minimal Left-Right symmetric extension of the SM. 
Remarkably, due to the relation between left and right Yukawa couplings in (3.6) -see also 
eqs. (3.8) and (3.9)- this triple vector correlation can be used to discriminate between 
charge-conjugation or parity as the Left-Right symmetry. More precisely, if the Dirac 
masses of neutrinos are smaller than the masses of the charge leptons, a value for the triple 
vector correlation bigger than 10“^ can only be the consequence of charge-conjugation as 
the Left-Right symmetry. 


A Kinematics of the /x —eq process and the triple vector correlation 


In this appendix we give some tools that could be useful when computing the triple vector 
correlation shown in Eq. (4.7) for the /x —>■ ey decay. 

For the anti-muon we use the spinor v{p^+) given by 


v{p^^+) 


-y/p^ i) ’ 


(AT) 


where = 1 and p^+ is given in Eq. (4.3). As shown in Fig. 1 the polarization vector of 
the muon is given by: 

s = |^(sin $ cos T, sin sin T, cos <1>) (A.2) 

and it is straightforward to show that in this case 


f e cos I \ 
y eW sin f J 


(A.3) 


One may find the same result by requiring ^ to be an eigenvector of a ■ h, where h is a 
unitary vector in the direction of s. 
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For the electron and for the reference frame shown in Fig.l we use 


Ve+{Pe+) = 


\Pe 


( —2e* 2“ sin ^ \ 
sin ^ 

• f) 

2 ie* 2 cos Y 
y —2e“*^ cos Y ) 


(A.4) 


The photon has two possible polarizations along the direction of motion and in the 
particular frame we are considering in Fig.l its polarization vector is given by, 




0 \ 

±i cos 
=Fi sin 05 
1 


(A.5) 


where we can explicitly see that when Qg = 0, the photon can only have a polarization 
±1 along the y-axis and and Pg+ are the 4-momentum of the outgoing photon and 
electron respectively - see Eq. (4.4) and (4.5). Once the expressions for the spinors of 
the participating fermions and the polarization vector of the photon are known, it is easy 
straightforward to compute the triple vector asymmetry given in (4.7). 

We found that the total decay rate is given by 


Ttotai — -|- |A/jp). (A.6) 

TT ^ 

It would be interesting to compare the above equation with the result one gets when 
summing the decay rates for cos > 0 to that of cos < 0, namely 

r(cos > 0) -|- r(cos < 0) = G^pm\e^{cos^ + sin^ (A.7) 

On the other hand, by subtracting the total decay rates for cos $ > 0 to that of cos $ < 0 
one gets: 

r(cos<I> > 0) — r(cos<I> < 0) = sm6s'Am{ALAp) (A.8) 

from which the asymmetry shown in (4.7) can be readily computed. 


B /i —)■ e total conversion rate and the triple vector correlation 
B.l Total conversion rate 

In this appendix we briefly comment about the amplitude of the /x —>■ e conversion process 
and the Born’s approximation we used. 

In computing the /x —)• e conversion process, one usually assumes the so called Born’s 
approximation for the outgoing electrons. This approximation has two meanings: one is 
computing the conversion rate to a given order in some small coupling; and the other is 
the assumption that electrons coming from the conversion process are plane waves. The 
point is that we can do better and have a complete control of both approximations at the 
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same time. More precisely for the relativistic one-electron atom and in the limit of big r 
(r >> ro, where V{r > rg) = 0), the solution of the Dirac’s equation at first order in the 
perturbation ffeff is of the form [16] 







V^TTxi^ip) 




(B.l) 


where 'ijji is any stationary state of the Coulomb potential, is one of the continuum 
energy solutions and H^ff is the effective Hamiltonian for the fi ^ e conversion process. 
Furthermore it can be shown that ipas is an eigenfunction of a ■ p + (3 with eigenvalue E 
so that ipas describes, indeed a plane wave [16]. In the high energy limit -neglecting the 
electron mass- the solution "ipas simplifies to 


pt'E' 

ipas = {ip'^\Heff\'ipi) 

K,fl 


Xk{p) \ 

-x^-.ip)) ■ 


(B.2) 


Finally if we are interested in computing the total conversion amplitude per unit flux (for 
a detector placed at fixed radius r = R) the total conversion rate is given by 


^conv 





2GU\Cl\^ + \Cr\^) (B.3) 


where the coefficients Cl and Cr are defined in section 5 and we may absorb the y/^ 
factor into the normalization of the wave function tppp in order to agree with the conventions 
adopted in [15]. 


B.2 Triple vector correlation in the conversion process 

In this appendix we give details of the calculation for the triplet correlation asymmetry 
in the /r —)• e conversion process within the formalism developed in [16]. Since we are 
interested in describing particles with a given polarization, we are going to make use of the 
spin projection operators for Dirac spinors. Instead of using the covariant spin projection 
operator we make use of the following projection operator 

pW = l(l±0-ho), (B.4) 

where 

O = p3a + {1 - p3){a ■ p)p (B.5) 

and no is the direction of the spin polarization vector in the rest frame of the particle, p 
is the direction of its momentum and the ± represent positive and negative polarization 
respectively. Furthermore it can be shown that the description of the spin with this operator 
is equivalent to the usual one given by the manifestly covariant spin operator Notice that 
the non-relativistic limit of can be taken in a transparent way by replacing /I —>■ 1. 

® see [17] chapter III. 
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For our present problem we assumed the muon to be non-relativistic and in the frame 
shown in Fig.l its polarization vector is of the form 

n^ = (0,ho), (B.6) 

where 

ho = (sin <i> cossin <i> sincos <i>) (B.7) 

by multiplying the wave function of the muon in the conversion process by one obtains 
the wave function of a non-relativistic muon with the given polarization. For the electron 
instead a full relativistic treatment is required since its energy is — e^, where 

is the muon mass and is the binding energy of the muon in the Is state of the muonic 
atom. In this case the spin projection operator coming from the conversion process is given 

by 


Pi+l = l(l + 0'.n5) (B,8) 

and 

O® • hg = /3(T • hg (1 - I3){a-pe){pe ■ h^), (B.9) 

ho = (0,1,0), pe = (sin6»s,cos6»s,0). (B.IO) 

Finally the wave function describing the polarized outgoing electron -coming from the 
conversion of a polarized muon- is obtained by applying to the solution (2.14) and a 
direct computation shows (for a detector placed at a fixed radius R): 

Wconv(cOs4' > 0) - a;conv(cos4> <0) = j dO • sgn(s^ • {pe X Se)) ■ 

= ^G|sind,3f?e[P(^-i-'5+B(CR - Cl){{C*r + Cl))] = G|sin0,9m(C'LC^) + 0{aZ) 

+ 0(^), (B.ll) 

since in the high energy limit the Coulomb phases satisfy 

<5-i-<^+i = | + 0(^). (B.12) 

The Coulomb phases (5±i are defined in Eq. (2.15) and dO. is given by dO. = dTd<l>sin4>. 
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